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Charge ordering in 0- (BEDT-TTF) 2 X materials 
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We investigate theoretically charge ordered states on the anisotropic triangular lattice char- 
acteristic of the (9-(BEDT-TTF)2X materials. Using exact diagonalization studies, we establish 
that the charge order (CO) pattern corresponds to a "horizontal" stripe structure, with ...1100... 
CO along the two directions with larger electron hopping (p-directions) , and ...1010... CO along 
the third direction (c-direction) . The CO is accompanied by co-operative bond dimerizations 
along all three directions in the highest spin state. In the lowest spin state bonds along the 
p-directions are tetramerized. Our theory explains the occurence of a charge-induced high 
temperature transition as well as a spin gap transition at lower temperature. 
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Fig. 1. Schematic structure of the 0-(ET)2X conducting layers, 
consisting of the ET cationic molecules. Solid lines correspond 
to stronger hopping tp, dashed lines to weaker hopping t c . 



Ground states involving spatially inhomogeneous 
charge distributions have in recent years been observed 
in a wide range of novel electronic solids, particularly in 
reduced dimensions. In the case of the organic charge- 
transfer solids (CTS), it has been suggested that the ap- 
pearance of charge-order (CO) may be related to super- 
conductivity. 1,2 ) In the present Letter, we focus on CO 
within one particular class of the CTS, the 0-(ET) 2 X ma- 
terials (here ET is short for BEDT-TTF). The 0-(ET) 2 X 
materials are particularly interesting for several reasons: 
(i) many detailed experimental studies have been per- 
formed; (ii) some members (e.g., the X = I3 salt 3 ' 4 )) of 
this class of CTS are superconducting; and (iii) the struc- 
ture of the ET conducting layer is fully two dimensional 
(2D), implying that standard single-particle theories of 
broken symmetry based on nesting are invalid. In this 
Letter we present exact numerical many-body calcula- 
tions establishing that the CO within the 6*-(ET) 2 X ma- 
terials is a combined bond-charge density wave (BCDW) 
state exhibiting a 2D extension of the "...1100..." type 
CO and BCDW that we have previously shown exist in 
a wide variety of quasi-one dimensional CTS. 1 - 1 The 2D 
bond distortions of this state also explain the spin gap 
seen at low temperatures in the f3-(ET) 2 X. 

The structure of the 2D conducting layer in the 9- 



type materials is shown in Fig. 1. We briefly review the 
experimental results, focusing on the most widely inves- 
tigated 9- (ET) 2 MM'(SCN) 4 materials. Here M is any 
of the monovalent cations Rb, Tl and Cs, while M' is 
one of the divalent cations Co and Zn. Slowly cooled 9- 
RbZn exhibits a metal-insulator (M-I) transition at 190 
K that is accompanied by a period doubling along the 
c-direction. 5 ' 6 ' Direct evidence for CO was found from 
NMR studies in this material, although the actual pat- 
tern of CO was not discussed. 7 ' From analysis of the 13 C- 
NMR line shapes, Chiba et al. have concluded that the 
CO in 0-RbZn corresponds to the horizontal stripe struc- 
ture (see Fig. 2(c)). 8 ' A similar conclusion was reached 
also from the analysis of the reflectance spectrum. 9 ' The 
190 K transition has no effect on the magnetic suscepti- 
bility, but a spin-singlet phase is obtained in this mate- 
rial below 30 K. 6 ' 7 ' 10 ' The lower temperature transition 
has been widely referred to as a spin-Peierls (SP) tran- 
sition. 6 ' 

In 0-RbCo, the M-I transition also occurs at 190 K and 
again there is a second transition (below about 40 K) to 
a spin singlet state. 5 ' 6 ' Crystal structure analysis gave 
clear evidence of dimerization along the c-direction 5 ' be- 
low this temperature. Importantly, since the dimeriza- 
tion in the c-direction has already occurred at the M-I 
transition this dimerization cannot itself be the origin of 
the spin gap. 

The 0-TlZn material is already semiconducting at 
room temperature, but a resistive anomaly is seen at 
165 K that may be accompanied by a lattice dimeriza- 
tion along the c-direction. 5 ' Evidence for short-range 
dimerization has been also found in 0-CsCo at the re- 
sistive transition, which, however, occurs at a much 
lower temperature of 20 K. 11 ' 12 ' In 6>-CsZn, NMR mea- 
surements have suggested a phase transition to a CO 
state near 20K. 13 ' Finally, a different 0-(ET) mate- 
rial, (9-(ET) 2 Cu 2 (CN)[N(CN) 2 ] 2 , is also of interest in the 
present context. This material shows a semiconductor to 
semiconductor transition at 220 K that is accompanied 
by a weak lattice dimerization along the c-direction. 14 ' 
No magnetic anomaly is seen at this temperature. Below 
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about 50 K there is a rapid drop in the spin susceptibil- 
ity, while the lattice dimerization along the c-direction 
gets stronger. 14 ) 

Summarizing the experimental data, the most notable 
features are (a) the existence of two transitions, one in- 
volving charge at high temperatures and the other in- 
volving spin at low temperatures; (b) where clear ev- 
idence exists, the horizontal stripe charge order (see 
Fig. 2(c)) is found; and (c) lattice dimerization along 
the c-direction occurs already at the M-I transition and 
therefore involves the charge degrees of freedom. We 
emphasize that the dimerization along the c-direction 
must necessarily involve distortions of intermolecular dis- 
tances, and not just charge modulation, since ordinary 
(as opposed to anomalous) X-ray scattering experiments 
determine only bond modulations and not site charge 
modulations. 15 ) The origin of the magnetic transition 
remains an open question, since the dimerization that 
appears at high temperature cannot itself produce a spin 
gap at lower temperature. 

Both Seo 16 ) and McKenzie et al. 17 ) have proposed 
models for the charge ordering in #-(ET). Following these 
authors and our own previous work, 1 ) we investigate the 
electronic properties of the ET layer using the extended 
Hubbard model: 

H = H Q + H ee (la) 
H Q = -J2 tijc{ a Cj, a (lb) 

H ee = Uy~] ni,im,i + V y~] mrij (lc) 

1 (ij) 

In the above, i,j are site indices, (...) implies nearest 
neighbors, nj = tij^ + nj t ±, and a is spin. The aver- 
age number of electrons (or holes) per site p equals \ 
for these ^-filled materials. In accordance with extended 
Hiickel calculations, 5 ) we consider t c < 0, t p > 0, and 
\t c \ « \t P \ within Eq. 1. However, we expect V p and 
V c to be roughly the same size, since while the hopping 
integrals depend upon the relative orientations of the 
molecules, the Coulomb integrals depend only on the dis- 
tances between them. Hence one must include both V p 
and V c . We have not included explicitly either intra- or 
inter-site electron-phonon couplings in Eq. 1. We will 
argue below that the signatures as well as the effects 
of such couplings can be anticipated from accurate so- 
lutions to the static extended Hubbard Hamiltonian for 
the ET lattice. 

Seo 16 ) has examined three possible CO patterns within 
the lattice of Fig. 1 at quarter-filling. All of these involve 
"stripes" , which can be vertical, diagonal, or horizontal 
stripes, as shown in Figs. 2 (a), (b), and (c). The lat- 
tice of Fig. 1 may also be thought of as a square lattice, 
with the essential proviso that hopping and Coulomb in- 
teractions along one of the diagonals must be included. 
In the square lattice picture, the bonds along the two 
p-directions occur along the x and y directions, respec- 
tively, and the intrastack bond along the c-direction in 
this case corresponds to one of the two diagonals of the 
square lattice. In this picture, the vertical stripe struc- 




Fig. 2. Different possible CO patterns in the 9-ET materials, (a) 
vertical stripe (b) diagonal stripe, and (c) horizontal stripe. 
Filled (unfilled) circles correspond to molecular site with greater 
(smaller) charge density. Corresponding stripe orders within a 
rotated rectangular lattice with hopping t p along the x and y- 
axes, and t c along the — x + y diagonal are shown in (d), (e), and 
(f). The hopping along the other diagonal is zero. 



ture becomes the simple checkerboard pattern of site oc- 
cupancy. Similarly, the diagonal stripe structure has al- 
ternate fully occupied and vacant rows along the x di- 
rection. Hence, both the vertical and diagonal stripe 
charge orderings are based on ...1010... CO in differ- 
ent lattice directions. In contrast, the horizontal stripe 
corresponds to a 2D version of ...1100... CO, 1 ' with 
the ...1100... orderings occurring along both x and y di- 
rections and ..1010... CO occurring along the diagonal 
correspond to the original c-axis. In nearly all cases that 
Seo studied using Hartree-Fock (HF)/Hartree solutions 
of Eq. 1, he found the vertical stripe to be the ground 
state, competing occasionally with the diagonal stripe. 
Only by introducing unequal intersite Coulomb correla- 
tions V p < V c and bond alternation along the c-direction 
did Seo find a CO pattern that is not vertical. In this 
specific case, the ground state was either the diagonal 
stripe or the horizontal stripe, but these were so close 
in energy that it was difficult to distinguish between the 
two cases (see Fig. 15 of reference 16 ^). McKenzie et al. 
studied the square lattice without t c or V c in the U — > oo 
limit within the slave-Boson approximation, 17 ) and for 
finite U with exact diagonalization. 18 ) Ignoring V c ob- 
viously strongly favors the checkerboard pattern CO in 
the square lattice (see Fig. 2(d)) which corresponds to 
the vertical stripe order (see Fig. 2(a)), and unsurpris- 
ingly these authors find vertical stripe CO in their model 
and assert that this is what is observed in the 6-ET CTS. 

Several important issues remain unresolved from these 
studies. Seo's work was done within the mean field ap- 
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proximation, 16 ' which is known to give erroneous behav- 
ior for the Hamiltonian of Eq. 1; for instance, in ID, the 
critical V for ...1010... CDW formation within HF theory 
is much too small, and the behavior of this transition as 
a function of U is incorrect. 19 ' 20 ' McKenzie et al.'s cal- 
culation 17 ' 18 ' ignores the Coulomb interaction V c com- 
pletely. While the hopping t c is indeed small and might 
legitimately be neglected, V is of similar magnitude in 
both p and c directions (because the distances are simi- 
lar), and V c /t c is actually much larger than V p /t p . Im- 
portantly, previous theories 16-18 ' have failed to explain 
(a) the c-axis dimerization (since the vertical stripe CO 
exhibits neither charge nor bond dimerization in the c 
direction), and (b) the low temperature spin gap transi- 
tion. It is argued in reference 16 ' that inclusion of e-ph 
interactions would lead to bond alternation among the 
occupied sites along individual ID stripes (i.e. along the 
...1-1-1... bonds in Fig. 2), as in a normal ID spin-Peicrls 
transition. Alternatively, Seo and Fukuyama have pro- 
posed a frustrated spin model to account for the spin 
gap. 21 ' However, it is not clear how the frustrated and 
CO states could coexist. 

To go beyond the previous approximations and as- 
sumptions, we have performed exact diagonalization 
studies for a 16-site lattice with 8 electrons, for the lat- 
tice structure of Fig. 1 (see Fig. 3 for the actual finite- 
size lattice that was investigated). We choose boundary 
conditions that are periodic along t c and both the t p di- 
rections, such that all three CO patterns of Fig. 2 arc 
well-defined and can be studied in an unbiased manner 
(i.e., the classical energies of the three stripe patterns 
obtained by setting tij — 0, are identical). We consider 
all hoppings t p to be identical, and similarly all hop- 
pings t c are taken to be the same. Since finite periodic 
systems do not exhibit broken symmetry, we follow our 
previously established procedure 1 ' and add a site energy 
component J^. to the Hamiltonian of Eq. (1), where 
ti is negative for the "occupied" sites and positive for the 
"unoccupied" sites. This is equivalent to including on- 
site electron-phonon (e-ph) coupling effects (within the 
classical approximation), with fixed spring constants for 
all three stripe structures. We then calculate the lowest 
energy corresponding to each specific CO for the small- 
est e, that still gives measurable energy differences be- 
tween the three stripe structures. As the spring-constant 
phonon energy is the same for on-site e-ph interactions 
for all stripe orders considered, the electronic energies in 
the static limit for small provide correct measures of 
the relative stabilities of the different CO patterns in the 
limit of + on-site e-ph coupling. 

We have chosen parameters similar to those in refer- 
ence, 16 ' and our results below are for U = 0.7 eV, t c — 
-0.01 eV and t p — 0.14 eV (note the overall minus sign in 
H in Eq. (1)), with V ranging from 0.15 - 0.35 eV and 
\ei\ = 0.01 eV. Our range of V covers the entire realistic 
region V < \U , given the fixed value of U = 0.7 eV. 
We have confirmed that increasing t c by a factor of 2 - 
3 or decreasing |ej| further does not lead to any modifi- 
cation of the phase diagram shown below. As we have 
discussed extensively in our previous work, 1 ' finite-size 
calculations tend to overemphasize the ...1010... CO and 




Fig. 3. The periodic 16 site lattice investigated for all three stripe 
CO orders. Charge densities and bond orders within the hori- 
zontal stripe order are shown for the ground spin singlet state 
for the parameters \t p \ = 0.14 eV, \t c \ = 0.01 eV, U = 0.7 eV, 
V = 0.35 eV, and |ej| = 0.01 eV. Note the bond dimerization 
along the c-direction and the bond tetramcrization along the 
p-directions. A similar calculation for the ferromagnetic state 
shows the same bond dimerization along the c-direction, but in- 
stead of tetramcrization, bonds are now dimerized even along 
the p-directions (see text). Differences in the bond orders oc- 
cur spontaneously as a result of charge ordering even though 
hopping integrals are uniform along t c and t p directions in our 
calculations. 



under- emphasize the ...1100... CO. Thus the relative en- 
ergy of the horizontal stripe CO in the present case is 
an upper limit, and the actual phase space region where 
this CO dominates in the thermodynamic limit is larger 
than what we find below. 

In Fig. 4 we have plotted the energy gained upon stripe 
formation, AE = E(a = 0) — E(a =/= 0) for the verti- 
cal, diagonal and horizontal stripe phases, respectively, 
versus V for the parameters given above. We find that 
in all cases the ground state contains either the diagonal 
stripe or the horizontal stripe. Importantly, the hori- 
zontal stripe CO dominates over a broad region of the 
phase space (V >0.18 eV), despite the isotropic intersite 
Coulomb integrals and uniform t c . 

In addition to the energies, we have calculated the 
charge densities rij and bond orders, the latter defined as 

the expectation value of bij = J2cr( c ia c j°- + c ] a c i^)- We 
find that the order parameter An here can be very large 
(~ 0.3), for |e,| = 0.01 eV, although this value will de- 
pend on the actual strength of the on-site e-ph coupling. 
However, the bond orders in a given direction are uni- 
form within both the vertical and diagonal stripe phases. 
On the other hand, within the horizontal CO, the bond 
orders alternate along the c-direction, while along the 
two p-directions there occurs a tetramerized periodic- 
ity. This is shown in Fig. 3. There is thus a strong 
tendency to spontaneous bond distortions in the hori- 
zontal CO phase in all three directions: the differences 
in the bond orders in Fig. 3 suggest that for arbitrar- 
ily small inter-site e-ph coupling there would occur bond 
dimerizations along the c-direction and bond tetramer- 
izations along the two diagonal p-directions in the hor- 
izontal stripe phase. The tendency to bond distortion 
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Fig. 4. The energy gained upon stripe formation, for the three 
stripe patterns of Fig. 2. Circles, squares and diamonds corre- 
spond to vertical, diagonal and horizontal stripe patterns, re- 
spectively. For V > 0.18 eV, the ground state has the horizontal 
stripe CO. 



in the 9-ET lattice is a consequence of the co-operative 
BCDW nature of the horizontal stripe phase. ^ If finite 
intra- and/or inter-site e-ph couplings are included, the 
cooperative nature of the BCDW accompanying horizon- 
tal stripe phase will further lower its energy. It is then 
likely that even in the small V region of Fig. 4 that the 
horizontal stripe will dominate. 

We have performed similar calculations for the ferro- 
magnetic spin configuration (total spin S = S max = 4), 
and have determined that the bond dimerization in the 
c-direction persists even here, although the magnitude 
of the dimerization is smaller. Thus, lattice dimerization 
in the c-direction is expected at the M-I transition, as 
soon as high spin states of the horizontal stripe order are 
formed. In contrast to the bond tetramerizations along 
the p-directions in the low spin state, however, the bonds 
in the p-directions in the high spin state are dimerized: 
the 1-1 and 0-0 bonds are the same, and the 0-1 and 1-0 
bonds are the same, but the two pairs are different. From 
our previous work, 1 ) there can be no spin gap in this S = 
Smax state with bond dimerizations, but spin gaps will 
appear in the S=0 "p-tetramerized" state along both p- 
directions. We further note that the horizontal stripes 
can not be viewed as ID half-filled systems, because the 
"unoccupied" sites in this phase actually have small spin 
components that are coupled to the large spin compo- 
nents on the occupied sites, rendering the spin behavior 
in this phase fundamentally 2D. 

In reference 16 ^ Seo found that explicit inclusion of a 
bond dimerization along the c-direction (in the form of 
alternating t c ) lowers the relative energy of the horizon- 
tal stripe phase (though not necessarily giving this as 
the ground state). The present exact calculations indi- 
cate that the lattice dimerization in the c-direction is a 
consequence of the horizontal stripe order, rather than 
being the cause of it. 

Our conclusions can be summarized as follows. We 



have shown that exact many-body calculations incorpo- 
rating all relevant electron-electron interactions predict 

horizontal stripe CO in the 6>-ET materials for a substan- 
tial region of parameter space. This result contradicts 
the previous results of Seo 16 ) and McKenzie et al 17 ' 18 ) 
and is consistent with experiments. In both the ferro- 
magnetic state and the low spin state the CO produces 
bond dimerization along the c-direction, with the ampli- 
tude of the dimerization larger in the low spin state. This 
explains the observed appearance of c-axis dimerization 
at the M-I transition 5,6 ' n ~ 14 ) as well as its increase at 
low temperatures, 14 ) where the free energy is dominated 
by the low spin state. In the two p-directions, bonds are 
dimerized in the ferromagnetic state but tetramerized in 
the low spin state. Our work thus explains naturally 
two distinct transitions in the #-ET materials, a high 
temperature transition 5, n ~ 14 ) that occurs as soon as the 
charge-occupancies correspond to the horizontal stripe 
CO, and a low temperature magnetic transition 6 ' 7,10 ) 
The spin gap transition is different from the usual SP 
transition, which is due to bond alternation in the ID 
half-filled spin system. The formation of a spin gap in 
2D in the absence of frustration is a unique feature of the 
strongly correlated 2D ^-filled band. Finally, our results 
raise questions regarding the theory of charge- fluctuation 
mediated superconductivity proposed by Merino et al., 2 ) 
as this work is based on the incorrect assumption of a 
checkerboard vertical stripe CO. 
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